Josephson-plasma and vortex modes in layered superconductors 
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The Josephson-plasma and vortex modes in layered superconductors have been studied theoretically 
for low magnetic fields parallel and perpendicular to layers. The two modes belong to the same 
lowest-frequency branch of collective-mode spectrum localized near vortices. This is the Josephson- 
plasma resonance if pancakes are strongly pinned and cannot move. Otherwise the lowest-frequency 
mode is a vortex mode governed by pancake pinning and the vortex mass. The latter is strongly 
enhanced by wandering of a vortex line around its average direction. The recently observed jump 
of the magnetoabsorption-resonance frequency at the vortex phase transition line is interpreted as 
a transition from the Josephson-plasma to the vortex mode. 

PACS numbers: 74.60.Ge, 72.30.+q, 74.72.Hs, 74.25.Nf 



Magnetoabsorption microwave resonances observed in 
Bi compounds (see and references therein) have be- 
come one of the most important subjects for studying 
the interlayer Josephson coupling. The resonances were 
observed for magnetic fields normal and nearly paral- 
lel to superconducting layers, and in the absence of the 
field. A remarkable feature of the resonances was that at 
high magnetic fields H the resonance frequency decreases 
roughly as 1/VH (the anticyclotronic behavior). Re- 
cently the magnetoabsorption resonances were observed 
in low perpendicular magnetic fields H . In low fields 
the resonance frequency very weakly depends on H, but 
at a field about a few hundred G, the magnetoabsorp- 
tion resonance frequency jumps to lower values with the 
anticyclotronic dependence, as observed in earlier exper- 
iments in higher fields. 

The first attempt to explain these resonances for high 
perpendicular magnetic fields was done by Kopnin et al. 
PJ who related the resonances to a vortex mode governed 
by surface pinning and the Magnus force. They were able 
to explain the anticyclotronic behavior, but the essential 
role of surface pinning and the Magnus force was not 
been confirmed by further experiments. 

Later the majority of researchers came to conclusion 
that the magnetoabsorption resonances are related to the 
Josephson Plasma Resonances (JPR) (see [pi and refer- 
ences therein) with the frequency 



(1) 



where u> p is the Josephson-plasma frequency at zero mag- 
netic field, f is the inplane coordinate, and fo{r) is the 
stationary gauge-invariant phase difference between lay- 
ers, which is nonzero due to misalignment of pancakes in 
neighboring layers by thermal fluctuations and disorder. 

However, this interpretation suffers from a number of 
inconsistencies (see discussion in Refs. [||-^|), and the 
original idea to relate the magnetoabsorption resonances 
with a vortex mode |I| should not be ruled out, though 
properties of this vortex mode must be different from 
those assumed in Ref. H. In Ref. || the magnetoab- 
sorption resonances in high parallel magnetic fields were 



interpreted as a vortex mode with the frequency deter- 
mined by the balance of the pinning force on pancakes 
and the inertia force determined by the vortex mass. The 
Magnus force is not effective in the parallel geometry. 

In the present work I investigate how a growing low 
magnetic field, cither parallel or perpendicular to layers, 
changes the collective-mode spectrum in layered super- 
conductors. The goal is to find correspondence of the 
JPR and the vortex mode to branches of this spectrum. 
In order to avoid a semantic ambiguity, I shall use the 
name JPR only for the mode with the frequency given by 
Eq. (Q), as done in all recent experimental papers, though 
in old papers this name had not the same meaning (see 
below) . The vortex mode is defined as a Goldstone mode 
which has a finite frequency only due to pinning. I shall 
show that the JPR mode and the vortex mode belong 
to the same branch of the spectrum which is the lowest- 
frequency oscillation mode localized near vortices. This 
mode becomes the JPR only if pancakes cannot move be- 
cause of very strong pinning. If the pinning strength is 
finite, the lowest-frequency mode is a vortex mode gov- 
erned by pinning and the vortex mass. The vortex mass is 
related with the electric energy in the interlayer spacing, 
and the misalignment of pancakes in the perpendicular 
field essentially increases the vortex mass. On the basis of 
the presented analysis I suggest that in the recent magne- 
toabsorption experiments in low perpendicular magnetic 
fields HH they observed a transition from the JPR to 
the vortex mode accompanied by a frequency jump. 

The analysis will be done for a single long Josephson 
junction, but results of the analysis are relevant also for 
a layered superconductor after slight modifications (see 
below). The state of the junction is determined from the 
sine-Gordon equation for the phase: 
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where c s = cy/ s/2e\ is the Swihart velocity, s is the 
barrier thickness, or the interlayer spacing, A is the Lon- 
don penetration depth for bulk superconductors forming 
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the junction (A 3> s), e is the dielectric permeability, 
Aj = (f>o/327r 3 A£ , j is the Josephson length, $ is the 
magnetic-flux quantum, and Ej is the Josephson cou- 
pling energy. 

The collective modes correspond to small oscillations 
around the stationary solution ipo(x): ip(x,t) — ipo(x) + 
ip'(x,t). The small phase ip' oc e~ luJt is determined by 
the equation (we shall omit the prime later on) 
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In the absence of the magnetic field the ground state 
corresponds to ipo — 0, and the solutions of Eq. (||) are 
plane waves exp(ifcr — iujt). The oscillation spectrum has 
a plasma edge, i.e., u> = ^Ju 2 + c s k 2 > lu p , where uj 2 = 

e 2 Ejs/eh 2 = c 2 /X 2 j is the Josephson plasma frequency. 

If the magnetic field is strictly parallel to the junction 
plane, the phase ip(x, t) depends only on one spatial coor- 
dinate x. In the stationary case (ip = 0) the sine-Gordon 
equation has exact periodic solutions ipo(x) which corre- 
spond to chains of Josephson vortices and the period 
(intervortex distance) a determines the magnetic field: 
H = $ /2Aa. For such tp (x), Eq. (§) can be also solved 
exactly for any intervortex distance ]Tc| ]. However, we 
introduce small periodic ipo(x) which is not a solution 
of the sine-Gordon equation and vanishes at distances 
more than Aj/2 from vortex centers. This yields a small 
periodic attractive potential U = [1 — cos po(x))/X 2 in 
Eq. (||). We assume that a 3> Aj (low magnetic fields), 
and look for a periodic solution symmetric in the interval 
—a/2 < x < a/2. The coordinate x = corresponds to 
the center of the potential well. At |x| < Aj/2 the per- 
turbation theory yields <p(x) oc 1 — J Q dx\ J 1 dx-iUix-i). 
At Aj/2 < \x\ < a/2 tp(x) oc cosfc(a/2 =p x). The lowest- 
frequency mode corresponds to the bound state with an 
imaginary k = ip. Continuity of the phase and its deriva- 
tive at x ~ ±Aj/2 requires that 

ptanhpa = —j / dx[l — cos <po(x)] . (4) 

* J J -a/2 

At very low magnetic fields when pa > 1, the bound 
state corresponds to the oscillation frequency 



Wo = <4 ~ c s P = u P \ 1 - [1 - (cos <po(x))] \ . (5) 



But if the inequality pa -C 1 holds, expansion in the left- 
hand side of Eq. (||) yields 

P « — 7j / ^L 1 ~ cosv? (a;)] = ^ — , 

aA J J -a/2 A J 

(6) 

and the squared mode frequency to 2 is given by Eq. ([!]). 



Thus the JPR belongs to the lowest-frequency mode 
localized near the potential well. However, the potential 
is very weak, the size of the bound state 1/p (the local- 
ization length) is larger than the intervortex distance a, 
and therefore localization is not so pronounced. 

Let us consider also the lowest-frequency oscillation 
which belongs to the continuum spectrum with a real k. 
The phase distribution for a continuum mode must have 
at least one node in the vicinity of the vortex, in order 
to satisfy the condition of orthogonality to the ground 
bound state. Thus roughly k — it /a, and the lowest fre- 
quency uj c of the continuum spectrum (delocalized phase 
oscillation) exceeds the zero-field plasma edge io p : 
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The frequency of this mode grows at increasing the mag- 
netic field. Note that in the old papers JTJ] this mode 
with lo c > lo p , but not the localized mode with loq < lo p , 
was called the plasma mode. 

Our calculation has not revealed the Goldstone vortex 
mode, because our potential does not correspond to a 
real Josephson vortex with the phase po which satisfies 
the stationary sine-Gordon equation. For real vortices 
the perturbation theory is invalid, but the exact theory 
yields that the lowest frequency luq vanishes. The phase 
in the bound state is proportional to the derivative of the 
stationary phase distribution d(po(x)/dx. 

So the procedure to tune the potential from zero to 
values of the order 1/Aj goes through states which are 
not realized physically, since the vortex is a topological 
nonlinear excitation, and its topological charge cannot 
change continuously. However, this procedure clearly 
demonstrates genesis of the JPR and the vortex mode 
from the original spectrum of the spatially uniform state. 
The JPR and the vortex mode belong to the same spec- 
trum branch related to a phase oscillation localized near 
the defect (vortex), but the JPR with the frequency Eq. 
([!]) is possible only for an artificial weak potential. Thus 
in parallel fields the JPR mode doesnot exist at all || . 

In a layered superconductor a perpendicular field cre- 
ates vortices which are chains of pancakes. If pancakes 
lie on an ideally straight line normal to layers, there is no 
phase difference across Josephson junctions. But because 
of disorder and thermal fluctuations, pancakes do not 
form an ideally straight line, and pancakes in neighboring 
layers are separated by random distances JHJ . For a sin- 
gle Josephson junction this means that the vortex line in 
two banks of the junction meets the junction plane in two 
different points with the distance r w between them, and 
a phase difference <fo(r) across the junction appears. If 
r w 3> A j, the phase distribution corresponds to a Joseph- 
son string, which is a segment of a Josephson vortex of 
width Aj stretched between the points (x = 0, y = 0) and 
(x = 0,y = r w ). If r w <C Aj, the statinary sine-Gordon 
equation yields that at r <C Aj <pn{r) is a dipole field: 
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At large distances r = \J x 2 + y 2 » r w (but not larger 
than Aj) the phase is <po(r) — r w x/r 2 . At distances 
r ^> A j the phase (fo(r) is exponentially small. We shall 
call such a phase distribution a short Josephson string. 

We look for the lowest-frequency mode for a short 
Josephson string assuming for the sake of simplicity that 
the potential —(1 — cos</3o)/A 2 in Eq. (|^) is axisymmet- 
ric. The most important perturbation originates from 
distances r w <C r <C A,/ where the axisymmetric part of 
1 — cost^oM ~ < Po( r ) 2 /2 is equal to r^/Ar 2 . An axisym- 
metric mode is determined from the equation 
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For a bound state in a weak potential one may neglect 
the right-hand side of Eq. (^) for r < A,/, and the per- 
turbation theory yields 
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Outside the potential well, at r > Aj, the axisymmet- 
ric bound-state oscillation is ip = AI (pr) + BK {pr), or, 
if pa <C 1, <p w A (l + p 2 r 2 /4) — Blnpr. The periodical 
boundary conditions in a real vortex array can be ap- 
proximately simulated by the condition <p'{r) — at the 
boundary of the Wigner-Zeitz cell r = a = \J <&q/ttH. 
Then 
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The continuity conditions at the border of the potential 
well (r — Xj) yield 
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In very low magnetic fields (a — > oo) 
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Then the squared frequency u 2 = lj 2 — c 2 p 2 is 
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For a long Josephson string r w 3> Aj one maynot use 
the perturbation theory, but the lowest-frequency mode 
can be calculated as a bending oscillation for an elastic 
string of the length r w and width ~ Aj with its ends 
fixed. The line tension of the string (the string energy 
per unit length) is e = 8Ej\j. The mass of the moving 
string is determined by the electric energy: 
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where E = (%/2es)<p — ~(h/2es)(vL-V)pQ is the electric 
field parallel to the c axis, and vl is the string velocity. 
Neglecting the logarithmic contribution from vicinities of 
the string ends, this expression yields for the mass per 
unit length fi t = M/r w = $§/2tt 3 c 2 sAj @. Finally, the 
frequency is u 2 = Tt 2 e/pir 2 w = it 2 u 2 p X 2 j / r 2 w = 7r 2 c 2 /r 2 . 
At r w ~ Aj where the crossover from a short to a long 
Josephson string takes place, this frequency is of the same 
order as given by Eq. (|l4|). 

We see that low magnetic fields do not affect the 
lowest-frequency localized mode. But at increasing H, 
i.e., decreasing the parameter pa in Eq. ([l2|), the loga- 
rithmic term in the denominator becomes unimportant. 
Then the squared frequency is given by Eq. (Tfl) , with 
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and the bound state corresponds to the JPR mode. The 
restriction on observation of the JPR mode is 
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Thus the JPR is a weakly localized mode with imaginary 
k, contrary to Bulaevskii et al. || who related the JPR 
to the most homogeneous delocalized state. Delocalized 
states belong to the continuum spectrum, and as well as 
in the parallel ID geometry [see Eq. ([?])], the lowest- 
frequency mode of the continuum spectrum exceeds the 
plasma frequency and grows with the magnetic field. 

Up to now we considered oscillations for a string with 
fixed (strongly pinned) ends, since if the ends moved, a 
singular contribution would appear in tp. Decreasing the 
pinning strength, the frequency of the localized mode 
decreases, and if pinning vanishes, the frequency must 
vanish also. The string ends correspond to pancakes in 
a layered superconductor. Thus the localized mode is a 
JPR if pinning of pancakes is infinitely strong. But for 
a finite pinning strength the localized mode is a vortex 
mode with pancakes and strings moving together and 
with the frequency which depends on pinning. Since 
pinning varies from a sample to a sample, it is difficult 
to suggest an universal estimation of the frequency of 
the vortex mode. But, as already discussed in Ref. M, 
the fact that the experimentally measured c-axis critical 
current has the same dependence on the magnetic field 
as the squared magnetoabsorption resonance frequency, 
confirms our suggestion that the magnetoabsorption res- 
onances in high fields correspond not to the JPR, but to 
the vortex mode governed by pinning. 

In order to observe a vortex-mode resonance in a lay- 
ered superconductor, the inertia force proportional to the 
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vortex mass (we neglect the Magnus force here) should 
not be small compared to the friction force. A moving 
vortex line is a chain of pancakes connected by random 
strings of average length r w , and r w may be consid- 
ered as a size of an extended vortex core. The Josephson 
length Xj = js (7 is the anisotropy ratio) usually ex- 
ceeds r w . The mass M of the short string is determined 
mostly by the vicinity of two pancakes, and in Eq. ( |l5| ) 
d(po/dx = —y/r 2 [or (y — r w )/r 2 ]. Then the mass per 
unit length is (cf. Ref. p[) 
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where £ a {, is the coherence length in the ab plane. On the 
other hand, the dissipation rate for string motion is also 
proportional to the electric energy: r\v 2 L = J o c E 2 (df)2, 
where a c is the normal conductivity along the c axis and 
77 is the friction coefficient per unit vortex length along 
the c axis. Then the quality factor of the resonance (the 
ratio of the inertia force to the friction force) is Q = 
jjLiOJo/r] — eujo/4:ir<7 c . It is identical to the quality factor of 
the uniform plasma oscillation as one can see comparing 
the supercurrent and the normal current in the expression 
for the total current j = —{eu> 2 /AiTiu})E + a c E. Thus 
normal-current dissipation is not more dangerous for the 
vortex mode than for the plasma mode. 

The assumption of infinitely strong pancake pinning 
cannot be valid everywhere, and therefore the JPR mode 
cannot be an universal explanation of the magnetoab- 
sorption resonances. This explanation was based on Eqs. 
(|l|) and (16). In Refs. ||,|8| I argued that in order to de- 
rive the observed anticyclotronic dependence from these 
equations, one must make a too strong assumption that 
r w ~ a in a wide interval of fields. The condition r w ~ a 
means a complete disintegration of vortex lines. Note 
that now Bulaevskii et al. || also believe that this con- 
dition does not take place and line structure retains in 
the vortex liquid. Then they must accept that it is im- 
possible to explain the anticyclotronic behavior in the 
liquid phase in terms of JPR. As discussed in Ref. [Q, 
the line structure in the vortex liquid was confirmed by 
experiments on the microwave absorption |l4f| . 

The analysis given above shows that Eq. ([j]) itself is 
not always valid. In addition, in high magnetic fields 
(a <C Aj) one may not use Eq. ( |l6| ) for the decoupled 
state of the vortex matter in which pancakes positions in 
neighboring layers are uncorrelated, but the 2D itralayer 
crystalline order still exists. Then the Josephson strings 
form an ordered 2D lattice in any interlayer spacing. In 
high magnetic fields when Aj exceeds a, Eq. ( |l6| ) is valid 
if orientations of the Josephson strings are random, and 
there is no interference between contributions of distant 
strings (at distances larger than a) to (ipo(r) 2 ) |ll|. But 
if the strings form a lattice, the distribution of fo{r) is 
similar to the distribution of the current component j y for 



the vortex lattice in the mixed state. This means that the 
contributions from distant strings to ipo(r) 2 cancel, and 
the logarithm \n(Xj/r w ) in Eq. ( |i"tj| ) must be replaced 
by \n{a/r w ). Then even the assumption r w ~ a does not 
help to obtain the anticyclotronic behavior. 

In summary, we explain the experimental dependence 
of magnetoabsorption resonances on a low perpendicular 
magnetic field |^,|| as the following. In very low fields, 
before the frequency jump, the magnetoabsorption res- 
onance weakly depends on a field as was predicted for 
the JPR ||, and therefore is the JPR mode. After the 
jump which apparently coincides with the vortex-matter 
phase transition (the second-peak line at low temper- 
atures), the magnetoabsorption resonance becomes the 
vortex mode governed by pancake pinning. This explains 
the anticycltronic dependence after the jump. The na- 
ture of the frequency jump depends on the structure of 
the vortex matter in the high-field phase which is not well 
established up to now. It maybe related with a jump of 
r w as supposed by the decoupling scenario of the phase 
transition [jl5| . The jump of r w should result in a jump 
of vortex mass, or in weakening of pinning, or both. 

Any scenario of the phase transition supposes that in 
the high-field phase directions of vortex lines strongly os- 
cillate jl6|, Then even the electric field strictly parallel 
to the c axis and to the average direction of vortices can 
make some tilted segments of vortex lines to oscillate. 
Thus the Lorentz-force-free geometry of the magnetoab- 
sorption experiments in perpendicular fields does not rule 
out a possibility to excite the vortex mode. 

The work was supported by the grant of the Israel 
Academy of Sciences and Humanities. 



[1 

[2 
[3] 
[4] 
[5 
[6 

ir 

[s 

[9 
[10 

[11 
[12 



[13] 
[14] 



M.B. Gaifullin et al, Phys. Rev. Lett. 83, 3928 (1999). 
T. Shibauchi et al, Phys. Rev. Lett. 83, 1010 (1999). 
Y. Matsuda et al, (unpublished). 
N.B. Kopnin et al, P hys. Rev. Lett. 74, 4527 (1995). 
L.N. Bulaevskii et al. 



:ond-ma t/9907462| . 



E.B. Sonin, Phys. Rev. Lett. 79, 3732 (1997). 
M. Gaifullin et al, Phys. Rev. Lett. 81, 3551 (1998). 
E.B. Sonin, Phys. Rev. Lett. 81, 3552 (1998). 
E.B. Sonin, Phys. Rev. B 60, 15 430 (1999). 
A. L. Fetter and M.J. Stephen, Phys. Rev. 168, 475 
(1968). 

A.E. Koshelev et al, Phys. Rev. B 53, 2786 (1996). 
The Josephson-string mass in low magnetic fields is field- 
independent, in contrast to the Josephson- vortex mass 
in high fields where the electric fields from neighboring 
vortices overlap {J. 

E.B. Sonin et al, Phys. Rev. B 57, 575 (1998). 

H. Enriquez et al, Phys. Rev. B 58, R14 745 (1998). 



[15] B. Horovitz, Phys. Rev. B 60, R9939 (1999). 
[16] P.H. Kes et al, J. Phys. I France, 6, 2327 (1996). 



4 



